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ABSTRACT 

A main disadvantage of using spectral methods for nonlinear conservation laws lies in 
the formation of Gibbs phenomenon, once spontaneous shock discontinuities appear in the 
solution. The global nature of spectral methods then pollutes the unstable Gibbs oscillations 
overall the computational domain, and the lack of entropy dissipation prevents convergences 
in these cases. 

In this paper, we discuss the Spectral Viscosity method, which is based on high frequency- 
dependent vanishing viscosity regularization of the classical spectral methods. We show that 
this method enforces the convergence of nonlinear spectral approximations without sacrificing 
their overall spectral accuracy. 
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1. INTRODUCTION 


We consider the 27r-periodic one-dimensional system of conservation laws 


(1.1a) 


du df{u) _ 
dt dx 


with prescribed initial conditions, u(x,t = 0) = Uo(x). We recall, [4], that solutions of (1.1a) 
may develop spontaneous jump discontinuities (shock waves) and hence the class of weak 
solutions must be admitted. Moreover, since there are many possible weak solutions, the 
system (1.1a) is augmented with an entropy condition which requires 


(1.1b) 


dU(u) dF{u) < 

dt dx 


Here, U(u) and F(u) = /“ < U'(w), f'(w)dw > is any entropy function and the correspond- 
ing entropy flux pair associated with (1.1a), so that a strict inequality in (1.1b) reflects the 
existence of physically relevant shock waves in the entropy solution of (1.1a), (1.1b). 

We want to solve the 27r-periodic initial-value problem, (1.1a), (1.1b), by (pseudo-) spec- 
tral methods. To this end, we use an N - trigonometric polynomial 


u N (x,t) = J2 u k (t)e ikx , 
k=-N 

in order to approximate the spectral or ^dospectral projection of the exact solution, u(x,t). 
Let Snu and 'ij)N u denote, respectively, the spectral-Fourier and the ^dospect rah Fourier 
projections of u(x,t), and let Ppju stands for either one of these two projections. Then, 
starting with u//(x,0) = P/yUo(x), the standard Fourier method, [1], lets evolves at 

a later time according to the (2 N + l)-dimensional approximate model 

( 1 . 2 ) = 

We can rewrite this approximation in the equivalent form 

(1.3) Eu n = (i- PnVM- 

The expression inside the right brackets is the local error (due to spatial discretization). The 
approximation (1.2) - or equivalently (1.3), is consistent with the conservation law (1.1) in 
the sense that its local error does not exceed 

(1.4) || Eu n = (I - PN)f(uN) || < Const-i • -^||uAr||tfi- 

In fact, the usual spectral estimates, [10], tell us that the approximation (1.2) is arbitrarily 
high-order (or spectrally) accurate approximation of the conservation law (1.1a). Namely, 
for any s > 1 there exists a constant, Const,, such that 


(1.5) 


\\Eu n \\ < Const, ■ N 'lluwll/r*- 
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The S V method was introduced in [1 1] where compensated compactness arguments were 
used to show convergence in the special case of the scalar Burgers’ equation. The conver- 
gence proof of the SV method was then extended in [5], [9], [12], to include general scalar and 
certain 2x2 systems of conservation laws. In the next sections, we outline the practical con- 
siderations involved in the implementation of the SV method, and we give a bird’s eye view 
of the above mentioned convergence results. Finally, we provide numerical evidence which 
shows how post-processing of the SV solutions enables to recover with spectral accuracy the 
pointwise values of the exact entropy solution. 

3. CONVERGENCE OF THE SV METHOD 

The modified equation associated with the SV method (2.2) takes the form 

(3-la) + ^[/(u„)] = 

The local error in this case, Eu consists of two contributions: the discretization error we 
had before in (1.3), and the additional viscosity error introduced on the right of (2.2), i.e., 

(3.16) Eun = (7 — P^)/(uat) + SnQn * 

Integrating (3.1) against u # over the 2x-period we obtain 
. 1 d r 2 * , , . , ,di I k ' ' 

(3- 2 ) j dtJo u N( x > t ) dx = H-far> Eu N)- 

Thus, the quantity on the right represents the amount of (quadratic) entropy dissipation 
rate. The counterexample discussed in Section 1 tells us that the control of such quantity is 
necessary for convergence. This brings us to 

DEFINITION. [12, Section 3]. The approximation (3.1a) is consistent with the entropy 
condition, if there exist constants, Const. > 0 and I 0 such that 

(3.3) -(^p Eun ^ ~ - Const - £ M\-^\\ 2 + IM| 2 , £n > jj. 

Using compensated compactness arguments, the main results of [12, Section 6] asserts 
that this kind of consistency together with L°°-stability imply convergence. We shall use 
this framework to prove the convergence of the SV method. 

The special form of the local error, Eun^ in (3.1b) reveals that the entropy consistency 
requirement (3.3) is fulfilled, if a sufficiently large amount of SV regularization, CnQn, is 
present. Indeed, taking into account the a^riori estimate (1.4) and using Parsevahs relation, 
we conclude that entropy consistency is achieved with SV kernels satisfying, 

(3.4) Q k (t) > Const. . 
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Can we use such a viscosity kernels without sacrificing spectral accuracy? using a vanishing 
viscosity amplitude of order, say, we find by (3.4), that the viscosity coefficients, 

Q k (t), should be activated only for high modes with wavenumbers \k\ > m N , where 


(3.5a) 





The resulting SV method then takes the form 

(3.56) + S~[PNf(uN)\ = jy H k 2 Qk{t)uk(t)e . 

It differs from the standard Fourier method (1.2) by the spectrally small amount of viscosity 
added on the right of (3.5b). Indeed, the contribution of this kind of viscosity to the local 
error in (3.1b), does not exceed 

||£v<?N * ^-|| < Const , • 

which retains the (formal) overall spectral accuracy of the Fourier method 

(3.6) \\Eu n \\ < Const a ■ 

Thus, by augmenting the Fourier method with a spectrally small amount of viscosity, we 
are able to enforce a sufficient amount of entropy dissipation so that the entropy consistency 
requirement (3.3) is met. Moreover, in the next section we show that this kind of spectral 
viscosity guarantees that our approximation, u N (x,t), remains uniformly bounded 

(3.7) !!«*(*, OIU- < M{t). 

Consequently, the (entropy) consistency and Testability imply that the SV method con- 
verges for a wide class of scalar conservation laws as well as certain 2x2 conservative systems, 

[9], [12]. 


4. DECAY AND CONVERGENCE RATES 

In order to gain a better insight into the role of spectral viscosity, it is instructive to 
study the decay rate of the Fourier coefficients, Uk(t) and fk{t) = //c(ujv(0)- This program 
was carried out in [5] for the special case of Burgers’ equation, and in [9] for the general 
scalar problem. An iterative argument outlined in [5], [9], shows that the presence of SV 
separates the computed spectrum into three different regions. 

1. Wavenumbers in the “inviscid” region \k\ < m^. The corresponding amplitudes, Ufc(f) 
and f k (t) are then governed by the underlying nonlinear conservation law. 

2. Wavenumbers outside the “inviscid” region, m N <\k\< N. Here we find that due to 
nonlinear interaction with the spectral viscosity, the corresponding amplitude dissipate 
at rate 

(4.1a) |«k(()| + l/i(<)l < Const, ■ (^p) rn„<\k\<N. 

In particular, we have here a transition to a third distinctive region, namely 
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3. Wavenumbers located at the highest portion of the spectrum, say < |jfe| < N. 
Then, by (4.1a), the corresponding amplitudes are negligibly small, i.e. , with vanishing 
viscosity amplitude of order e N >N~ 2 &, we have 


(4.16) |tijb(f)| + |A(<)I < Const, • N (1 20)t + e~ Nt , P <\, \n < |Jfe| < N. 

We conclude that the presence of spectral viscosity enforces a spectral decay of the 
discretization error 

(4.2) ||(7 - Pjv)/(wjv)|| < Const, • + e~ Nt , P <\- 

We observe that the spectral decay of the discretization error is valid independently whether 
the underlying solution is smooth or not, [5, Section 3], [9, Section 3]. This shows that the 
SV solution is essentially governed by the equation 

(4 .3) 

This equation is closely related to the standard viscous regularization (2.1). We note, how- 
ever, that unlike (2.1), the viscosity regularization in (4.3) is nonlocal due to the finite 
support of the convoluted kernel Qn(x,i). 

One can use now the viscous equation (4.3) in order to obtain a uniform bound on the 
SV solution. Integrated against pu p I ^ 1 (x, t) over the 2n- period, equation (4.3) yields 

C 4 - 4 ) fc\\ u N\\Lr < £n\\^{I - Qn]j^\\lp' 

Standard trigonometric estimates can be used to upper bound the RHS of (4.4), 

(4-5) e*||^[J " ||lp < Const.\\u N \\„. 

Indeed, the inequality (4.5) is just an Zr p -version of the entropy consistency estimate stated 
in (3.3), consult [12, Section 5]. Hence, by combining (4.4), (4.5) and carefully iterating on 
the I^-norms of us{x,t), we derive the 7°°-bound of uji/(x,t) promised earlier. 

Furthermore, Schochet [9] used the viscous equation (4.3) in order to conclude an almost 
optimal L x -convergence rate of the SV method, namely 


(4-6) \\u N - u\\ L i < Const.N-P, P < i 

z 

Figures 4.1 and 4.2 compare the behavior of the spectral-Fourier method with and without 
spectral viscosity in the case of Burgers’ equation, which is subject to initial condition, 
u 0 (x) = sinx. The resulting ODE system for the Fourier coefficients 


(4.7) 


^“*(0 + \ ik X) «p(t)u,(<) = -ek 2 Q k u k (t), 

p+q=k 


1*1 < N, 
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was integrated up to time t = 1.5, using the fourth-order Runge-Kutta method. The method 
was complemented with viscosity coefficients Q k = 1 only for wavenumbers \k\> where 
uiff ~ 2 N* . Figure 4.1 shows that the SV solution convergence strongly (but not uniformly) 
to the exact entropy solution, in sharp contrast to the oscillatory behavior of the viscosity- 
free Fourier method in Figure 4.2. Improved results were obtained in Figure 4.3, by using 
C°° viscosity coefficients, Qk, connecting wavenumbers in the inviscid region, \k\ < m N , 
and the highest wavenumbers, ] A:| ~ N. This kind of smoothly varying SV prevents the 
propagation of the Gibbs phenomenon into the whole computational domain that can be 
noticed in Figure 4.2, consult [11]. Moreover, Figure 4.4 shows this SV solution after it 
was post-processed by the spectrally accurate smoothing procedure discussed in [2]. Finally, 
in Table 4.5, we quote from [6], the pointwise errors of the post-processed SV solution for 
the Burgers’ equation. The results indicate the spectral convergence rate obtained by the 
post-processed SV solution in the shock-free zones of the entropy solution. 
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Figure 4.1. The SV method for Burgers’ Figure 4.2. The spectral Fourier method 

equation. for Burgers’ equation. 
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Figure 4.3 


Figure 4.4 


The pseudo-SV method with smooth viscosity kernel before and after post-processing. 
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|«at( 

x,t) - 

-u(x,t)\ 

at t = 1.5 

X — 

N-- 

= 16 

N 

= 32 

.314 

5.7- 

10" 3 

2.3 

• lb- 4 

.628 

3.9- 

lO" 3 

2.8 

• 10~ 4 

.942 

7.2- 

10- 3 

2.7 

• lO" 4 

1.25 

1.2- 

10- 2 

2.8 

l 

o 

T““ f 

1.57 

1.4 • 

To - 4 

2.8 

10~ 4 

1.88 

1.7- 

IF 5- 

3.0 

10" 4 

2.19 

2.0- 

IF 5- 

3.6 

10" 4 

2.51 

2.1 • 

hF” 

4.3 

•10~ 4 


Table 4.5 

Pointwise errors of the post-processed pseudo-SV approximation compared with the exact 
entropy solution of Burgers’ equation. 
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